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1 .  INTRODUCTION 


This  first  annual  report  describes  work  completed  by  Yap 
Analytics,  Inc.  under  contract  F19628-89-C-0033  during  the  period 

t 

from  April  20,  1989  through  April  30,  1990.  Our  objective  is  the 
study  of  the  mechanisms  involved  in  the  emissions  of  the  infrared 
(IR)  radiation  of  the  upper  atmosphere  and  the  related  impacts  of 
this  background  clutter  on  surveillance  sensors.  During  this  time 
frame,  we  concentrated  our  efforts  on  the  investigation  of  the  mo¬ 
dels  for  non-eguil ibrium  phenomena  leading  to  emission  from  very 
high  rotational  levels.  One  critical  parameter  in  the  model  is 
the  scattering  cross  section  in  a  complex  multi-body  collision  and 
interaction.  In  this  investigation,  the  impulse  approach  (IA)  is 
being  applied  to  atom-molecule  collisions  which  accounts  for  the 
total  scattering  amplitude  of  an  atom  incident  upon  a  complex  sys¬ 
tem  by  summation  of  independent  two  particle  interactions.  A  re¬ 
sultant  momentum  distribution  of  each  constituent  particle  is  then 
evaluated.  Studies  to  date  by  various  investigators  uses  simpli¬ 
fications,  the  peaking  approximation  (PA) ,  to  evaluate  the  complex 
integral  for  the  momentum  of  non-interaction  (spectator)  parti¬ 
cles.  Our  work  investigated  the  validity  and  accuracy  of  the  PA. 

During  this  period,  we  also  worked  on  atmospheric  code  vali¬ 
dation  and  compared  code  predictions  with  field  measured  data.  In 
particular,  we  compared  outputs  of  the  SHARC  code  with  the  RAD 
code  to  understand  the  models  used  in  each  code,  the  discrepancies 
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of  the  models  and  the  resultant  differences  in  the  code  predic¬ 
tions.  Adjustments  of  universal  constants  were  made  to  the  codes 
for  comparison  on  an  equal  basis.  A  standard  atmospheric  profile 
was  used  as  input  to  SHARC  code  to  provide  atmospheric  radiance 
profile  predictions  for  the  NO,  03,  and  C02  bands.  Comparisons 
of  these  predictions  were  then  made  with  field  data  measured  on 
the  SPIRIT  and  SPIRE  programs. 


2 


2.  ATMOSPHERIC  MODELING  AND  ANALYSIS 

The  Impulse  Approach  (IA)  to  collisions  was  proposed  as  early 
as  1959.  The  theory  suggests  the  total  scattering  for  a  projec¬ 
tile  incident  upon  a  complex  system  is  the  sum  of  the  two-body  am¬ 
plitude  between  the  projectile  and  each  of  the  constituents  of  the 
complex  system.  Thus  far  we  have  confined  our  work  to  atom-diatom 
collisions.  Previous  computations  have  always  been  accomplished 
by  using  an  additional  approximation,  called  the  peaking  approxi¬ 
mation  (PA),  to  greatly  simplify  the  formalism.  We  have,  however, 
developed  an  exact  impulse  approach  that  does  not  rely  on  this  ap¬ 
proximation.  The  results  obtained  using  our  exact  IA  programs 
have  been  used  for  comparison  with  PA  computations  and  with  vari¬ 
ous  experimental  data.  Additionally,  the  validity  of  IA  has  been 
determined  by  the  properties  of  semi-detailed  balancing  and  the 
optical  theorem.  These  findings  are  detailed  in  the  two  papers 
included  in  Appendix  A  and  B  of  this  report.  Impulse  formalism  for 
atom-molecule  collisions:  Inadequacy  of  the  peaking  approximation, 
published  in  the  August  1,  1989  issue  of  Physical  Review  A,  and 
Criteria  for  applicability  of  the  impulse  approach  to  collisions, 
which  will  appear  in  the  next  issue  of  the  same  publication. 

2.1  Comparison  of  exact  and  PA  results 

Differential  and  total  cross  sections  have  been  computed  for 
various  systems  to  evaluate  the  validity  of  the  peaking  approxi- 
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Taylor  expanded  t-matrix  and  the  peaking  t-matrix.  Figure  (3) 
superimposes  the  cross  section  ratios  of  Figure  (1)  with  the  ratio 
of  the  squares  of  absolute  values  of  the  analytical  t-matrix  and 
the  peaking  t-matrix.  It  is  clear  that  the  analytical  expression 
for  exact/PA  provides  an  excellent  approximation  over  nearly  the 
entire  range. 

We  have  further  examined  the  validity  of  this  analytical 
treatment  by  comparing  the  exact  and  differential  cross  sections 
with  those  obtained  using  the  analytical  approach.  Figure  (4) 
displays  the  differential  cross  section  based  on  the  exact  and 
analytical  approach  for  the  Li*  +  N2  process  (0,0)  -*  (1,6)  at  1 
eV  as  a  function  of  scattering  angle.  The  two  curves  are  nearly 
indistinguishable  except  for  around  14°  where  the  cross  section 
almost  vanishes. 

In  addition,  we  have  tested  the  validity  of  the  analytical 
approximation  for  many  other  systems,  including  H  +  N2  and  Li'  + 

N:  using  various  initial  and  final  molecular  states. 

2.2  Assessing  the  validity  of  IA 

Using  the  semi-detailed  balance  (sdb)  condition,  we  have  as¬ 
sessed  the  criteria  for  the  applicability  of  the  impulse  approach. 
The  sdb  condition  requires  the  cross  sections  in  the  forward  di¬ 
rection  and  the  normalized  cross  sections  in  the  reverse  direc¬ 
tions  to  be  equal.  A  violation  of  the  sdb  condition  indicates  an 
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internal  inconsistency  in  IA,  leading  us  to  believe  the  approach 
is  not  valid.  We  have  computed  cross  sections  usinq  our  exact  im¬ 
pulse  approach  for  various  systems  in  the  forward  and  reverse 
directions,  in  particular  Li*  +  N2  ♦-*  Li*  +  N,  (v',j').  Using 
the  half-on-the-shell  (post  and  prior)  models  of  the  two  body  t- 
matrix,  we  have  shown  that  in  both  the  IA  leads  to  a  violation  of 
the  sdb  conditions  for  small  scattering  angles.  The  cross  sec¬ 
tions  are  studied  quantitatively  as  a  function  of  the  relative 
translational  energy  and  the  mass  of  the  incident  particle.  For 
large  scattering  angles  we  also  find  that  the  cross  sections  are 
independent  of  the  choice  of  energy  parameter  6S  (the  kinetic 
energy  of  the  colliding  particles  before  the  collision) .  We  sug¬ 
gest  that  two  criteria  be  used  for  assessing  the  validity  of  IA: 

1)  the  sdb  condition  and  2)  model  independence,  where  each  model 
is  distinguished  by  the  choice  of  es  in  the  two  body  t-matrix. 

The  findings  are  detailed  in  our  paper,  Criteria  for  applicability 
of  the  impulse  approach  to  collisions  (Appendix  B) . 

In  addition  to  the  sdb  condition,  we  have  also  used  the  opti¬ 
cal  theorem  to  determine  the  validity  of  our  current  approach. 

This  involves  computing  total  cross  sections  (integrated  over 
scattering  angle)  for  all  possible  transitions  from  a  given  ini¬ 
tial  state  and  a  fixed  relative  translational  energy.  This  has 
been  done  for  both  the  Li’  +  N2  and  H  +  N2  systems  using  (v=0,  j=0) 
or  (v=0,  j=12)  as  the  initial  state  and  1  eV,  2  eV,  and  4  eV  as 
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the  incident  energies.  The  sum  total  of  cross  sections  is  then 


compared  to  the  elastic  scattering  amplitude  in  the  forward  direc¬ 
tion.  Our  findings  suggest  that  the  multiple  collision  terms  may 
be  necessary  to  satisfy  the  sdb  and  optical  theorem  conditions. 
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Figure  1.  Differential  cross  section  ratio  for  1  eV  and  various  scattering  angles 
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Figure  3.  Theoretic.nl  fit  of  the  squares  of  the  absolute 

the  differential  cross  section  ratio. 


as  the  only  emitter.  The  vibrational  temperatures  for  the  626 
isotope  v2  (01101)  level  were  compared  with  the  output  from  the 
RAD  program.  At  the  lower  altitudes,  60  km  to  100  km,  the  differ¬ 
ences  between  the  two  programs  were  relatively  small  and  the  va¬ 
lues  oscillated,  sometimes  SHARC  had  higher  values  but  more  often 
they  were  lower.  At  altitudes  higher  than  100  km,  however,  the 
differences  became  increasingly  larger  with  the  values  obtained  by 
SHARC  being  the  lower  of  the  two.  At  the  highest  altitude,  175 
km,  SHARC  was  24.6  K  lower  than  the  vibrational  temperature  deter¬ 
mined  by  the  6-1  code  (RAD) .  The  rate  coefficients  for  the  C02 
transition  00001  -»  01101  were  scrutinized.  The  expression  found 
in  SHARC  for  the  N2  and  02  molecules  rate  constant  was  similar 
to  the  expressions  found  in  RAD.  However,  the  rate  constant  in 
SHARC  for  the  O  molecule  was 

lxlO'^T172  +  2 . 32xlO~V76  75/rU3  +  7.41xl0''V82Q/T  (1) 

while  RAD  used 

1. 5xlO'l3Tl/2  +  2 . 32x10'V767  5/t1/3  (2) 

Notice  that  for  the  dominant  term  of  the  expression  SHARC  had 
lxio13  whereas  RAD  had  1.5xl0'13  as  the  coefficient.  This 
coefficient  in  the  expression  for  the  rate  constant  in  SHARC  was 
changed  to  match  that  of  RAD  and  SHARC  was  rerun.  The  same  trend 
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was  seen  as  in  the  first  run  but  the  temperatures  were  in  closer 
agreement.  Again  the  SHARC  temperatures  were  lower  and  the  most 
significant  difference  occurred  at  the  highest  altitude.  At  175 
km  the  temperature  difference  was  14.3  K. 

Due  to  the  continuing  discrepancy  between  the  codes  the 
Einstein  A  coefficients  were  compared.  New  coefficients  for  SHARC 
were  generated  and  the  C02  linking  file  was  changed  to  replace 
the  coefficients  with  the  new  ones.  The  new  Einstein  A  coeffi¬ 
cients  were  compared  with  those  found  in  the  open  literature  and 
also  with  those  in  RAD.  Most  values  compared  favorably  and  only 
one  Einstein  A  coefficient  was  found  to  effect  this  band  emission 
and,  therefore,  it  was  examined  more  closely  .  The  values  were 
close  enough  and  should  not  cause  such  a  large  discrepency  theo¬ 
retically  between  the  vibrational  temperatures  but  SHARC  was  run 
again  for  the  same  case  with  the  new  Einstein  A  coefficients. 

These  output  values  were  very  similar  to  the  previous  run.  The 
largest  discrepancy  occurred  again  at  the  highest  altitude  where 
SHARC  was  14.7  K  lower  than  RAD. 

In  order  to  compare  the  excitation  rates  SHARC  was  run  for 
an  upwelling  earthshine  case.  SHARC  does  not  include  in  its  out¬ 
put  the  collisional  excitation  rates  but  it  does  give  the  radia¬ 
tive  excitation  rates  for  the  earth,  sun,  and  atmosphere.  RAD's 
output  is  given  as  a  combination  of  these  radiative  excitation 
rates.  Therefore,  for  a  meaningful  comparison  the  sum  of  the  in- 
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dividual  SHARC  outputs  was  compared  with  RAD's  output.  Although 
discernable  pattern  emerged,  many  of  SHARC' s  values  were  about 
half  of  those  determined  by  the  RAD  code. 

The  three  C02  input  files  used  by  SHARC  were  limited  to  on¬ 
ly  the  information  necessary  for  this  one  reaction.  The  Einstein 
A  pertinent  to  this  reaction  was  set  to  1.525  for  both  codes.  The 
reaction  rate  in  SHARC  for  N2/02  is 

lxlO'l6T1/2  +  6 . 81xlO‘8e‘82/T'/3  (3) 

and  Peter  Wintersteiner  of  ARCON  rearranged  the  expressions  in  the 
RAD  code  to  match  this.  The  first  coefficient  of  1.5E-13  was  re¬ 
tained  for  the  O  for  reaction  rate.  The  problem  still  remained 
and  the  difference  in  temperature  at  175  km  was  11.0  K  with  SHARC 
being  the  lower  of  the  two.  Figure  5  shows  these  results. 

One  reason  for  the  discrepencies  in  the  vibrational  tempera¬ 
tures  is  that  SHARC  uses  a  Doppler  lineshape  and  RAD  uses  a  Voigt 
lineshape.  Emission  that  occurs  at  the  lower  altitudes  can  be  ab¬ 
sorbed  at  the  higher  altitudes.  The  Voigt  lineshape  has  wings  to 
account  for  this  but  the  Doppler  lineshape  does  not.  Since  the 
Doppler  lineshape  does  not  give  the  proper  emission  at  the  lower 
altitudes  the  absorption  at  the  higher  altitudes  is  too  small.  The 
Doppler  shape,  therefore,  is  giving  lower  temperatures  than  the 
Voigt  shape. 

The  radiative  pumping  for  each  layer  as  seen  from  different 
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observation  points  in  the  atmosphere  is  given  in  RAD's  output. 
Looking  at  each  observation  point  it  can  be  determined  which  la¬ 
yers  are  producing  the  most  pumping  at  each  height.  At  the  lower 
altitudes,  the  two  layers  closest  to  the  altitude  produce  the  most 
pumping  but,  at  higher  altitudes  the  largest  contribution  is  from 
the  total  pumping  below  the  lowest  altitude  (the  blackbody  contri¬ 
bution)  . 

Figure  6  depicts  the  radiative  pumping  rate  determined  by  RAD 
at  150  km  tangent  height  as  a  function  of  altitude  in  fractional 
form.  Each  rate  is  divided  by  the  total  absorption  at  this  height 
which  is  1.93  x  104  photons/  ( sec-cm3)  .  The  total  combined  con¬ 
tribution  from  below  40  km  as  a  fraction  of  the  total  is  labelled 
as  the  blackbody  contribution  in  Figure  6  and  is  the  largest  con- 
tributer  to  the  total  pumping  at  150  km.  The  blackbody  contribu¬ 
tion  is  about  5.8  %  of  the  total  contribution.  At  high  altitudes 
there  is  more  absorption  than  at  altitudes  that  fall  in  the  middle 
(80-100  km)  range.  These  middle  range  altitudes  have  a  lower  tem¬ 
perature  which  causes  a  narrow  absorption  cross-section  and  the 
blackbody  contribution  at  these  heights  would  be  much  less  than  at 
150  km.  The  Voigt  lineshape  allows  for  this  blackbody  contribu¬ 
tion,  however,  it  is  also  important  to  note  that  while  these  radi¬ 
ative  pumping  values  obtained  from  RAD  give  a  very  good  indication 
of  the  importance  of  the  wider  wings  it  is  still  not  an  ideal 
model . 
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Radiative  Pumping  Rate  at  150  km 
as  a  function  of  altitude 
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Figure  6.  Radiative  punning  rate  predicted  by  RAD  code. 


(lu>1)  1H0I3H  1N33NV1 


2U 


25 


RADIANCE  (W/cm2  -sr) 

Figure  10.  Comparison  of  SHARC  prediction  with  SPIRE  and  SPIRIT  data  for 
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Impulse  Formalism  For  Atom-Molecule  Collisions;  Inadequacy  of  the 
Peaking  Approximation. 

Ramesh  D.Sharma3,  Pradip  M.Bakshib,  and  Joseph  M.Sindonic. 

a.  Optical  Physics  Division,  Geophysics  Laboratory, 

Hanscom  Air  Force  Base,  MA.  01890. 

b.  Physics  Department,  Boston  College,  Chestnut  Hill,  MA. 02167. 

c.  Yap  Analytics  Inc.,  594  Marret  Road,  Lexington,  MA.  02173. 

ABSTRACT 

Expressions  for  differential  and  total  cross-section 
for  atom-diatom  scattering  are  derived  using  the  Impulse 
Formalism  without  any  approximations.  Results  for  the 
rotational-vibrational  scattering  are  obtained,  for  the  first 
time,  without  using  the  peaking  approximation  (PA) .  For  the 
specific  case  of  a  hard  core  potential,  it  is  shown  that,  except 
for  elastic  scattering,  PA  results  are  substantially  different 
from  the  true  impulse  results. 
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The  basic  idea  of  Impulse  Approach  (IA)  to  collisions  was 
proposed  by  Chew1.  The  total  scattering  amplitude  for  a 
projectile  incident  upon  a  complex  system  is  taken,  in  IA,  to 
be  the  sum  of  the  various  two-body  amplitudes- (the  projectile 
plus  each  of  the  constituent  particles  of  the  complex  system) ; 
the  constituents  not  involved  in  the  two-body  scattering  remain 
unaffected  during  the  impact  and  are  termed  spectators.  The 
role  of  the  binding  potential  of  the  complex  system  is  to 
generate  a  momentum  distribution  of  the  constituent  particles. 

IA  has  been  applied  to  atom-diatom  collisions  by  Bogan2, 

Eckelt,  Korsch  and  Philipp3'6,  and  by  Beard  and  Micha7. 

In  all  of  these  papers  an  additional  simplification,  the 
Peaking  Approximation  (PA) ,  is  invoked  to  evaluate  the  integral 
over  the  spectator  momentum.  In  this  approximation,  the  t- 
matrix  representing  the  two-body  scattering  process  is 
evaluated  for  a  particular  value  of  the  spectator  momentum. 

This  value  of  the  spectator  momentum  is  taken  to  be  one  for 
which  the  product  of  the  initial  and  final  state  molecular 
wavefunctions  is  maximum.  Since  the  two  body  t-matrix  varies 
much  more  slowly  with  internal  molecular  momentum  than  does  the 
product  of  the  molecular  wavefunctions,  the  integral  over  the 
spectator  momentum  factors  into  two  parts  :  (i)  two  body  t- 
matrix  evaluated  for  a  specific  value  of  the  spectator 
momentum,  and,  (ii)  the  molecular  form  factor  integral  which 
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can  be  evaluated  easily  for  a  given  value  of  the  momentum 
transferred.  PA  not  only  simplifies  the  computations  but  it 
also  gives  a  very  simple  physical  picture  of  the  impulse 
formulation.  However,  its  validity  and  accuracy  have  not  been 
quantitatively  assessed. 

We  have  recently  developed  a  formalism8  which  permits 
computation  of  IA  cross-sections  without  resorting  to  PA.  This 
allows  us  to  investigate  the  validity  and  assess  the  accuracy 
of  PA.  The  purpose  of  this  Rapid  Communication  is  to  point  out 
that  PA  is  valid  only  in  the  limited  range  of  vibrational ly  and 
rotationally  elastic  scattering.  It  gives  grossly  erroneous 
overestimates  for  vibrationally  inelastic  forward  scattering. 

It  also  does  not  converge  to  the  true  IA  results  even  for  high 
relative  translational  energies. 

The  state-to-state  differential  scattering  cross- 
section  for  molecule  1-2  to  undergo  transition  from  initial 
vibration-rotation  state  vj  upon  collision  with  atom  3 
to  final  state  v'j'  is  given  by35,8 


aS(vjp3-v'j'p'„9)=  - rT2jViT 


It 


where 
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(2) 


<$'3|  Tls)!  <J>3>=  J  dq3  <J>''  (q'3)<q's  I  tls)  |  q$>0  (  q3) 
and  |  <*>3>  p  |v,j,m,p>  is  the  initial  state,  the  final 
state  being  denoted  by  primes.  p3  and  p'3  are  the 

momenta  of  the  incident  particle  before  and  after  the  collision 
in  the  center-of-mass  (CM)  frame  of  the  atom-molecule  system. 

P3  is  also  the  momentum  of  the  incident  particle  3  with 

respect  to  the  CM  of  the  molecule  1-2.  p^  denotes  the 
momentum  of  particle  a  with  respect  to  the  CM  of  be  and 
q  is  the  relative  momentum  of  b  and  c.  This  set  of 

a 

momenta  are  called  Jacobi-momenta3,5 .  ^denotes  the 

reduced  mass  of  the  system  (a, be) .  0  is  the  scattering  angle, 

i.e.,  the  angle  between  p3  and  p3'.  Summation  over 

m  and  m'  has  removed  the  dependence  of  the  differential  cross- 

section  on  the  azimuthal  angle.  T(sl  and  t(sl  are  respectively 

the  three  body  and  two  body  transition  matrix  elements,  s  being 

the  spectator  atom.  $(q3)  and  (q3)  are  the 

initial  and  final  state  molecular  wavefunctions  in  the  momentum 
representation . 

We  have  shown8  how  the  integral  in  Eq.  (2)  can  be 
evaluated  without  any  approximations.  We  start  by  expanding  the 
two  body  t-matrix  element  in  a  spherical  harmonics  expansion  in 
the  unit  vector  q3, 
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with  t'5^  given  by 


tlSLM  (q3f*  )  =  Jd<33Y*LM  (0(3)  <q's|  tlsl|qs>.  (4) 

Substituting  Eq.(3)  in  Eq.(2),  after  some  algebra,  we  get 

«p  '  |T(sl|<f>  >=  I  (qi)°i,J'M  Y*Cy(q)  11%^  (q) 

pxypLM  3 

*C ( j  LX ; 00 ) C ( j  LK ; mMp. )  C ( \(3  j  •  ;00)C(Xpj  •  ;pym'  )  ,  (5) 

where  q  means  +  sign  for  s=l  and  -  sign  for  s=2,  q  is  the 
momentum  transferred  during  the  c  ^uion,  jp(2j+l),  etc., 

C's  are  the  Clebsch-Goruan  coefficients9,  and 

N(s)LMxp(q)=Jn"dr^2jp(  I «s  1  rq)  K'Wr)  Xv.r(r)  .  (6) 

•  S 

Wlth  ^p(-}  (i^)-  and 

K'SU(^)=  l/Jdq,  q^  jx(q3r)  t(slw(q3)  Ivj(q3),  (7) 

-VJ  (<q3)  =  r2  xvj(r)  j*j(q3r)f  (8) 

Xvj(r)  being  the  radial  part  of  the  vibration-rotation 
wavefunction  and  j  is  the  spherical  Bessel  function  of  order 


To  simplify  the  remaining  algebra  it  is  convenient  to 


choose  the  z-axis  of  the  coordinate  system  along  q,  the 
direction  of  momentum  transfer.  Then 


(9) 


Using  the  properties  of  Clebsch-Gordan  coefficients3,  C's,  and 
Racah  coefficients,  W's,  Eq. (5)  can  be  rewritten  as 

«t>  •  |T(S,|*>=  £  (qi)P  PNl:JtMxj}(q)  C(  j  LK;  00) 

J  J  pxu  4n j  '  p 

XC ( K(J j  '  ;00)C(  jJj  '  ;mM)C(L(U;M0)W(  jLj  *p;KJ)  .  (10) 

The  C  and  W  coefficients  express  the  relation  between  the  total 
angular  momentum  change  during  the  collision  J  and  the  initial 
and  final  rotational  quantum  numbers,  j  and  j ' .  We  also  note 
the  total  angular  momentum  change  J  is  composed  of  two  parts,  L 
coming  from  the  two  body  t-matrix  and  P,  the  usual  component, 
derived  from  the  momentum  transferred  during  the  collision. 

The  expression  for  the  differential  cross-section  is  now 
obtained  by  summing  the  absolute  square  of  the  collision 
amplitude  over  m  and  M=m'-m  leading  to, 


♦v '  j 


’  P  '  3  »  0)  =  (¥)* 


p' 

3 

*7 


WJM') 

I 

J=lj-J'l 


s 

M 


If;mI 


(11) 


where , 
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F^=  S(i)U"x*W  P  [|^]2  C(jLX;00)C(\pj  •  ;  00)  C  ( L()J  ;M0 ) 

xW(jLj'D;KJ)  [N(V(i+  (-D^Ul  (12) 

In  contrast,  the  differential  cross-section  in  the  peaking 
approximation  is  given  by3 


[ig(vjp3-vj'p,),e)]p,.  (^V 

x  |  fC»  +-11)  +f_ii4J  ft2)  ft2) 

I  1  V'J',VJ;X  ^  PA  ^  X  '  L  VJ-,VJ;XU  PA 


2 


V  \C2(jkj';00) 

X-.lj-j'l 


(13) 


where 


f  vr.vj;xPJ0dr  r2xvT(r)  jx(  |<*s|  rq)XVJ(r) 


(14) 


The  total  cross-section  is  obtained  from  the 


differential  cross-section  by  the  relation 


o(vjp-v'j'p')=^4  (d^)qdq- 


(15) 


Eq.(ll)  is  our  result  for  an  impulse  calculation 
without  further  approximations,  whereas  eq.(13)  is  the  standard 
PA  result.  To  compare  the  two  sets  of  calculations  we  pick  a 
hard  core  potential  to  represent  the  two-body  interaction.  The 


two-body  t-matrix  for  this  potential  is  available  in  a  closed 
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form10  and  is,  in  fact,  the  only  potential  used  in  the 
previous  atom-diatom  studies.  We  compare  the  results  for  a 
system  extensively  studied  previously3,6 

Li*+N2(v,j)  -  Li*+N2(v',j').  (16) 

The  computations  have  been  performed  for  initial  vibrational 
level  v-0  and  final  levels  v'=0,l,  and  2.  For  relative  kinetic 
energy  of  1  eV,  the  hard  core  radius  is  taken  to  be3 
1.62A.  In  PA,  for  homonuclear  diatomic  molecules,  the 
product  of  the  initial  and  final  wavef unctions  is  assumed  to 
peak  at  the  mid-point  of  their  two  centers,  q3=0  and  q3'=0, 

and  the  two-body  t-matrix  is  evaluated  at  fixed  internal  momentum 
q3=tq,  where  q=  P3'-p3 

is  the  momentum  transferred  during  the  collision. 

Molecular  wavefunctions  used  in  this  calculation  were 
obtained  by  numerical  integration  of  the  radial 
Schrodinger  equation11  using  RKR  diatom  potential 
constructed  from  spectroscopic  constants12 . Previous  work 
used  harmonic  oscillator  potential  to  obtain  the  diatom 
wavefunctions.  There  are  substantial  differences  between 
peaking  calculations  using  harmonic  oscillator  and  more 
accurate  wavefunctions  for  situations  involving  large 
vibrational-rotational  inelasticity.  To  highlight  the 
differences  between  the  peaking  and  exact  IA  calculations  we 
used  more  accurate  wavefunctions  throughout  this  work. 
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Fig.l  gives  a  plot  of  the  differential  cross-section 
for  the  process  (0,0-*l,6)  as  function  of  scattering  angle  0. 

For  forward  scattering  the  exact  IA  result  is  about  two  orders 
of  magnitude  smaller  than  the  PA  result!  Even  though  the  peaks 
and  valleys  for  the  two  calculations  track  each  other  closely 
the  discrppancy  between  the  two  sets  of  calculations  is  visible 
even  on  the  log  scale.  The  quantitative  difference  between  the 
two  calculations  is  best  described  by  the  ratio  p  of  the 
exact/PA  results,  displayed  in  fig.  2.  It  is  clear  from  this 
figure,  which  also  displays  this  ratio  for  (00-*0j'),  j'=0,10, 
and  20,  that,  except  for  elastic  (00-*00)  scattering,  there  are 
noticeable  differences  between  the  two  calculations.  These 
differences  appear  to  increase  with  energy  loss,  independent  of 
whether  the  energy  ends  up  in  vibration  or  rotation.  Further, 
even  for  large  angles  PA  overestimates  the  exact  results  by 
about  15%. 

Fig.  3  displays  the  total  cross-section 
aT(00-*n'j')  for  n'=0,l,  and  2  as  function  of  energy  loss. 

PA  grossly  overestimates  the  true  result  for  small  i '  for  n'=l 
and  2.  This  is  because  small  Aj  scattering  is  essentially 
forward  scattering  for  which  the  two  calculations  give  very 
different  results.  The  total  cross-section  mirrors  this 
disparity,  weight  factor  sin  0  notwithstanding.  For  large  j  * , 
when  almost  all  the  energy  loss  goes  into  molecular  rotation 
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EXACT  vs.  PEAKING 


SCATTERING  ANGLE 
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ABSTRACT 

Using  an  exact  formulation  of  impulse  approach  (IA)  to  atom- 
diatom  collisions  we  assess  its  internal  consistency.  By  com¬ 
paring  the  cross-sections  in  the  forward  and  reverse  directions 
for  the  vibrational-rotational  inelastic  processes,  using  the 
hal f-on-the-shel 1  (post  and  prior)  models  of  the  two-body  t- 
matrix,  we  show  that  in  both  cases  the  IA  leads  to  a  violation  of 
the  semi-detailed  balance  (sdb)  condition  for  small  scattering 
angles.  An  off-shell  model  for  the  two-body  t-matrix,  which  pre¬ 
serves  sdb,  is  shown  to  have  other  serious  shortcomings.  The 
cross-sections  are  studied  quantitatively  as  a  function  of  the 
relative  translational  energy  and  the  mass  of  the  incident 
particle,  and  criteria  discussed  for  the  applicability  of  IA. 

PACS  numbers : 03 . 80 . +r ,  34.50.-s. 
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Introduction:  Since  its  inception  the  Impulse  Approach  (IA) 
to  collisions  has  been  used  in  diverse  fields  of  Physics1-5.  The 
conditions  of  its  validity  are  often  heuristically  described, 
but  quantitative  measures  of  its  applicability  are  not  available. 
IA  as  applied  to  atom-molecule  collisions  to  date4'6  consists  of 
three-steps:  (1)  The  multiple  collision  expansion  is  truncated, 
retaining  only  the  single  collision  terms.  (2)  The  three-body  T- 
matrix  element  is  replaced  by  a  two-body  t-matrix  element.  The 
choice  of  the  energy  parameter  of  the  two-body  t-matrix  distin¬ 
guishes  the  various  models  (e.g.  post4,  prior4,  full5)  .  (3)  A 

further  approximation,  the  Peaking  Approximation  (PA) ,  was  in¬ 
voked  in  the  earlier  studies4,5  to  evaluate  the  transition 
ampl itude . 

We  have  recently  shown6,7  how  to  evaluate  the  transition  am¬ 
plitude  without  any  approximations,  eliminating  the  errors  due  to 
the  last  step.  This  now  allows  us  to  take  a  critical  look  at  the 
choice  of  the  energy  parameter  in  the  two-body  t-matrix.  In  this 
paper  we  show  that  the  two  hal f-on-the-shell  models4,  prior  and 
post,  violate  the  semi-detailed  balance  (sdb)  condition  for  small 
scattering  angles,  indicating  an  internal  inconsistency  in  those 
approaches.  We  also  show  that  the  (off-shell)  full  model5,  which 
satisfies  sdb,  has  other  serious  shortcomings.  We  further  suggest 
alternate  models  which  satisfy  sdb  (and  do  not  have  other  formal 
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difficulties)  to  test  the  model  dependence  of  the  cross-sections. 
Finally,  we  discuss  criteria  for  the  applicability  of  IA.  The 
detailed  results  given  here  are  obtained  using  a  two-body  inter¬ 
action  with  a  hard  core  potential1'5'8.  The  conclusions  reached, 
however,  should  be  applicable  to  other  potentials  as  well. 

In  a  multiple  collision  expansion4  of  the  atom-diatom 
three-body  T-matrix, 

T(z)  =  Tm+  T(2i+  T(iiG3T(21+  T(21G3T(i!+  ...,  (1) 

where  G3=  (z-H3)"'p  (z-  H0-  V3) -l  is  the  propagator  corresponding 
to  the  atom-diatom  Hamiltonian  which  includes  the  relative  trans¬ 
lational  energy  j  potential  energy  V3  of  the  molecule  1-2  but 
not  the  inter  ~tion  energy  of  the  incident  atom  with  atoms  1  (V  ) 
and  2  (V,  of  the  diatom.  The  first  two  terms  correspond  to  the 
collision  of  the  incident  atom  3  with  atoms  2  and  1  respectively, 
while  atoms  1  and  2  are  the  spectators,  indicated  by  super¬ 
scripts.  The  next  two  terms  represent  the  double  collision  terms. 

IA  retains  only  the  two  single  collision  terms1'6.  The 
second,  and  crucial,  step  in  IA  is  the  reduction  of  the  three- 
body  T-matrix  to  the  two-body  t-matrix  by  assuming  that  the 
time  duration  of  the  collision  is  much  shorter  than  the  cha¬ 
racteristic  times  for  the  molecular  motion.  The  function  of  the 
intramolecular  potential  is  then  to  generate  a  momentum  distri¬ 
bution  for  the  two  atoms  constituting  the  diatom.  Formally,  the 
two  steps  are  represented  by  ’ 5 
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«t>'3\  T  |  <t>  3>=>  «t>'3\  Tll)+  T(2)  |  <f>>, 


(1) 


and 

<0’3I  Tls)(E)  |  <t>  3>=>  f  dq3  0'’(q'3)<q's  |  t(s)  <  €s)  |qs>$(q3)  /  (2) 

where  |  <f>3>  p  |v,j,m,p>  is  the  initial  state,  the  final  state 
being  denoted  by  primes.  E,  the  eigenvalue  of  H3,  is  the  total 
energy  of  the  initial  (and  final)  state;  es=Es=  q2/(2mlu) 
is  the  relative  kinetic  energy  of  the  non-spectator  t  and  u 
atoms  before  the  collision  ,  p3  and  p'3  are  the  momenta 
of  the  incident  particle  before  and  after  the  collision  in  the 
center-of-mass  (CM)  frame  of  the  atom-molecule  system.  T(sl 
and  tlsl  respectively  generate  the  three-body  and  two-body  tran¬ 
sition  matrix  elements,  s  being  the  spectator  atom.  <*>(q3)  and 

$'( q3)  are  the  initial  and  final  state  molecular  wave- 
functions  in  the  momentum  representation.  q$  and  q 

are  the  Jacobi  momenta19  before  and  after  the  collision. 

The  choice  €S=ES  (the  kinetic  energy  of  the  colliding  pair 
before  collision)  in  the  two-body  t-matix  element 

<qs'  I  t(sl(€s)|qs>  =  t(s)(q^'  ,qs;€s)  (3) 

seems  to  be  the  natural  one  in  the  spirit  of  IA  and  is  called 
the  post  form4  of  IA.  Other  choices  have  been  considered;  £s  = 
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E's  =  q'2s/(2mlu)  is  called  the  prior  form'1,  m,u  being  the 
reduced  mass  of  the  t  and  u  atoms.  The  choice5  €s  =E  -  p‘/2ps, 
being  the  reduced  mass  of  the  s  and  (t+u)  system,  may  be 
called  the  full  form  in  that  it  subtracts  the  spectator  atom 
energy  from  E,  the  total  physical  energy  of  the  process.  The 
post  and  prior  forms  are  hal f-on-the-shell ,  since  one  of  the 
momenta,  q^  or  q  •;  lies  on  the  energy  shell.  The  full  form 

is  off-shell.  The  basic  assumption  of  IA,  that  the  spectator 
momentum  p^  has  not  changed  during  the  collision  is  incorpo¬ 
rated  in  each  case.  Calculation  of  cross-sections  using  different 
€.  will  generally  give  different  results.  We  return  to  this 
point  later. 

To  compare  the  relative  merits  of  different  models,  we  apply 
symmetry  considerations.  Since  the  Hamiltonian  is  independent 
of  time  and  quadratic  in  momenta,  it  is  invariant  under  time  re¬ 
versal.  Combined  with  the  space  inversion  symmetry,  this  leads  to 
the  semi  detailed  balance  (sdb)  relation10,11  between  the  dif¬ 
ferential  cross-sections  for  the  forward  and  reverse  processes, 

p[do(i-»f)/dfi]/p(f)=  p' [do(f-*i)/dfl]/p(i)  (4) 

where  p ( i )  and  p(f)  are  the  densities  of  state  in  the  initial 
and  final  states;  p  and  p'  represent  the  incident  flux  densities. 
A  valid  theory  must  lead  to  cross-sections  which  satisfy  Eq.  (5). 

Results ;  We  have  computed  the  differential  cross-sections  for 
the  much  studied  process'*  b 
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(5) 


L. *+  N2(  v, j)  ~  Li*+  N2(v' , j  ') 

in  the  forward  and  reverse  directions  using  IA,  Eq.  (3),  with 
the  previously  employed*1'5  hard  core  two-body  (atom-atom) 
potential12  for  various  values  of  the  relative  translational 
energies  Etel.  In  Figure  1  the  differential  cross-section  using 
the  post  model,  for  the  process  v=0  j  =  12  -*  v=l  j=6,  in  the 
forward  direction  at  Erel=  lev  is  given  by  the  solid  line.  The 
cross-section  for  the  reverse  process  multiplied  by  p  '  3p  ( f ) /  — 
p.p(i),  the  ratio  of  the  densities  of  state  times  the  inci¬ 
dent  flux  (p(f)a  [ ( 2 j '  +  1 ) p  '  )  and  p(i)a  ((2j  +  l)p]),  and 
termed  'normalized'  reverse  differential  cross-section,  is  given 
by  the  dotted  line.  The  two  curves  should  be  identical  if  the 
theory  obeys  sdb.  Up  to  a  scattering  angle  of  about  50°  the  two 
curves  do  not  appear  related.  From  50°  to  70°  the  two  curves 
have  a  similar  structure.  For  larger  scattering  angles  they  are 
indistinguishable . 

The  dotted  line  in  figure  2  gives  a  plot  of  the  ratio  of 
the  two  differential  cross-sections  of  Figure  1  as  a  function  of 
the  scattering  angle  for  Erf|  of  1  ev,  its  departure  from  unity 
providing  a  quantitative  measure  of  the  violation  of  sdb.  The 
solid  line  represents  a  similar  ratio  for  Erfl=  4  ev.  It  is 
clearly  seen  that  the  two  curves  which  appear  unrelated  in 
Figure  1  show  a  definite  periodicity  in  their  ratio,  and  the 
whole  structure  gets  squeezed  towards  smaller  angles  at  higher 
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relative  translational  energy. 

We  conclude  that  the  post  form  does  not  satisfy  the  sdb 
condition  for  small  angles.  If  we  use  the  prior  form  for  e., 
the  forward  process  is  represented  by  the  dotted  line  in  Fig.  1 
and  the  normalized  reverse  process  is  represented  by  the  solid 
line.  Thus  the  prior  form  also  fails  to  satisfy  sdb  where  the 
post  form  fails.  It  can  be  shown  from  the  structure  of  the  post 
and  prior  t-matrix  elements  that  the  post*-*prior  interchange 
is  equivalent  to  the  forward*-*  (normalized)  reverse  interchange. 

The  full  form  obeys  sdb,  a  result  predictable  from  the  structure 
of  its  t-matrix. 

Now  we  examine  the  mass  dependence  of  the  sdb  violation. 
Figure  3  shews  the  ratio  of  the  forward  and  normalized  reverse 
differential  cross-sections  for  the  same  N;  transition  caused  by 
a  collision  with  H  atom  (solid  curve)  or  Li*  (dotted  curve) , 
both  at  Erel  of  1  ev,  as  functions  of  r.q,  q  =  |  p  1 -p  |  being 

the  momentum  transfer.  This  is  a  more  convenient  parameter 
to  compare  processes  with  different  projectile  masses.  For 
r,q>  20  the  two  curves  show  a  period  of  n  with  the  H  atom 
curve  having  a  smaller  amplitude.  The  effect  of  a  lighter  inci¬ 
dent  particle  appears  similar  to  that  of  increasing  the  energy  of 
relative  motion  as  seen  in  Figure  2.  Figure  4  displays  (rcq),, 
the  smallest  value  of  rcq  at  which  sdb  is  satisfied  to  1  %  as 
function  of  relative  energyi:  for  the  two  incident  particles 
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studied. 


Discussion;  We  have  shown  above  that  both  the  post  and  prior 
forms  violate  the  sdb  criterion  for  small  angles.  The  other  model 
used  in  the  literature5,  the  full  form,  does  not  violate  sdb 
but  poses  other  problems.  The  first  is  a  conceptual  diffi¬ 
culty;  the  choice  es  =  E  -  p2^/2 always  leads  to  a  negative 

energy  €s  for  a  part  of  the  range  of  integration  over  q^  in 
Eq . ( 3 ) .  Thus  the  amplitude  integral  in  Eq.(3)  remains  ill- 
defined  without  introducing  a  cut-off  (in  the  range  of  integra¬ 
tion)  or  attempting  analytical  continuation  into  negative  ener¬ 
gies.  Secondly,  even  if  a  cut-off  is  introduced  as  a  practical 
measure,  the  full  form  leads  to  very  high  cross-sections  for  the 
smaller  angles,  yielding  an  unreasonably  large  total  cross- 
section  in  significant  excess  of  the  quantum-mechanical  shadow. 
The  total  cross-sections  based  on  the  post  and  prior  forms,  on 
the  other  hand,  understate  the  shadow  limit. 

In  view  of  these  intrinsic  difficulties  with  the  models  used 
in  the  literature,  one  can  consider  other  ad  hoc  choices  for  €s 
which  avoid  some  of  these  problems.  An  obvious  choice  is  £s  = 

^(Es  +  E's),  which  may  be  termed  the  mean  model.  This 
introduces  a  symmetry  in  the  structure  of  the  t-matrix  and  the 
sdb  criterion  is  automatically  satisfied.  It  is  also  free  from 
any  definition  problems.  Figure  5  provides  a  comparison  of 
these  four  models  for  the  same  process  as  in  Figure  1  and  2  with 
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Er„,  =  4  eV.  The  full  model  is  evaluated  by  cutting  off  the 


range  of  where  es  would  become  negative.  For  this  energy, 

the  wave  functions'  product  in  Eq. (3)  in  the  excluded  domain  is 
quite  insignificant,  making  this  a  practical  calculation.  (For 
an  energy  such  as  Erel  =  1  eV,  the  cut-off  for  the  full  model 
would  take  away  a  significant  part  of  the  form  factor  integral, 
making  it  a  meaningless  calculation.)  It  is  clear  from  Figure  5 
that  for  larger  angles  these  results  are  model  independent. 

For  intermediate  angles,  the  post,  prior,  and  mean  models  trace 
similar  patterns.  For  small  angles  the  models  vary  significant¬ 
ly- 

In  conclusion,  we  suggest  (1)  the  sdb  condition  and  (2) 
model  independence  as  two  criteria  to  be  employed  in  assessing 
the  applicability  of  IA.  In  the  domain  I  where  both  the 
criteria  are  met,  one  may  use  IA  with  confidence.  Figure  4 
provides  the  boundary  of  this  domain  for  1%  accuracy  for  both 
the  post  and  prior  models.  In  the  domain  II,  where  the  post, 
prior,  and  the  mean  models  retain  the  same  patterns  with  slight 
quantitative  differences,  IA  may  be  expected  to  be  a  good 
approximation  with  any  of  these  models,  with  possible  errors  of 
the  same  magnitude  as  their  differences.  The  full  model  seems 
to  overstate  the  cross-sections  and  seems  much  less  reliable  in 
this  domain.  For  the  small  angles,  domain  III,  all  the  models 
give  rather  different  results.  This  is  also  the  domain  where 
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FIGURE  CAPTIONS 


FIG . 1 .  Direct  (solid  line)  and  normalized  reverse  (dotted 
line)  differential  cross-sections  for  Li*  +  N2(v=0,j=12)  -»  Li* 

+  N2  (v '  =1 ,  j  '  =6)  at  relative  translational  energy  Ere,=  1  eV  as  a 
function  of  scattering  angle  0°. 

FIG. 2.  Ratio  of  the  forward  and  normalized  reverse 
differential  cross-sections  of  Fig.  (1)  as  function  of 
scattering  angle  9°  at  Efel  =  leV  (dotted  line)  and  Erel  =  4  eV 
(solid  line) . 

FIG. 3.  Ratios  of  the  forward  and  normalized  Reverse 
differential  cross-sections  for  the  (0,12)  <=>  (1,6)  transitions 
of  N:  at  relative  translational  energy  of  1  eV  for  collision 

with  Li*  (dotted  line)  and  H  (solid  line) . 

FIG.  4.  Threshold  value  of  rcq  for  1%  deviation  from  semi- 
detailed  balance  condition  as  function  of  relative  translational 
energy  in  eV  for  the  (0,12)  <=>  (1,6)  transitions  of  N:.  Dotted 
line  is  for  collision  with  Li*  and  solid  line  is  for  collision 
with  H. 
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FIG. 5.  Comparison  of  the  differential  cross-sections  for 
the  process  Li'+  N2(v=0,j=12)  -*  Li*  +  N;(v  '  =  1 ,  j  '  =6)  at  relative 
translational  energy  of  4  eV  for  four  models  of  £s: 
post(solid  line),  prior(dashed  line),  mean (dash-dot-dash)  and 
full (tiny  dashes). 
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FIGURE  B- 
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FIGURE  B- 


